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The corrections to g-mode frequencies caused by the presence of a central magnetic field and 
rotation of the Sun are calculated. The calculations are carried out in the simple one dimensional 
magnetohydrodynamical model using the approximations which allow one to find the purely analyt- 
ical spectra of magneto-gravity waves beyond the scope of the JWKB approximation and avoid in 
a small background magnetic field the appearance of the cusp resonance which locks a wave within 
the radiative zone. These analytic results are compared with the satellite observations of the g-mode 
. frequency shifts which are of the order one per cent as given in the GOLF experiment at the SoHO 

' board. The main contribution turns out to be the magnetic frequency shift in the strong magnetic 

field which obeys the used approximations. In particular, the fixed magnetic field strength 700 KG 
results in the mentioned value of the frequency shift for the g-mode of the radial order n=-10. The 
' rotational shift due to the Goriolis force appears to be small and does not exceed a fracton of per 

cent, an < 0.003. 



I ^! I. INTRODUCTION 

I The recent progress in the search for low frequency solar oscillations - possibly g-modes - achieved in the GOLF 

O (SoHO) experiment [1] and future experimental proposals [3] open new possibilities for the detailed studies of the 

-j-J . solar radiative zone (RZ) as well as the solar core which cannot be done by solving inverse helioseismic problem for 

^ ' the high frequency p-modes. 

I— I i The frequencies found in the experiment [l| turn out to be slightly shifted 6uj/uj ^ 1 % compared to the predictions 

• of the Solar Seismic Model (SSeM) 0, i] obtained in the absence of magnetic field and solar rotation. 

" ^' [ The influence of the RZ magnetic field on the g-modes frequency shifts was studied in Ref. 0] . In that work the 

I calculations were performed in the simplified one dimensional (1-D) model that allowed one to evaluate the scope of 

] the magnetic field perturbation theory without g-modes absorption (and reflection) by Alfven or cusp resonances layers 

^\ inside the RZ. This simpliflcation guarantees that g-modes can enter the solar surface. Such 1-D calculations appeared 

CS| to be in a good agreement with the later magnetic shifts estimates in three dimensional (3-D) MHD model [g. 

ly-^ ^ The correct separation of the cusp (Alfven) resonances in 1-D model (see Sec. IIII 5] below), which enables one to 
simplify the problem in question, is the signiflcant advantage in comparison with 3-D model that allows us to compute 

[ — , rotation and magnetic corrections. Such corrections are calculated in a 3-D model with help of the perturbation theory, 

> without referring to the presence of MHD resonances, using eigenfunctions of horizontal and vertical displacements 

J> ! computed in SSeM not accounting for a magnetic field or a stellar rotation. On the contrary, the eigenfunctions 

k> ] of the 1-D problem, the vertical velocities ^^"'''^'^"^ = —iu!{n,l,fl,Bo)^i^^'^'^''^"\z), depend on both the magnetic 
d ' and the stellar rotation. These functions are the exact analytical solutions to the considered problem. 

' In the present work we continue to study the modifications of analytical MHD spectra for internal helioseismological 
waves in the 1-D model beyond the JWKB approximation taking into account both the magnetic field and the rotation 
of the Sun. Note that in the 3-D MHD the solar rotation can be taken into account in analytical spectra only in the 
frame of the JWKB approximation of short waves, L <C 0^/5 f^. 

It should be also mentioned that in the 1-D model we loose the effect of the stellar rotation with respect to an 
observer because of the simplification of the topology. In the 3-D model this effect is characterized by the azimuthal 
frequency splitting ~ nifl. However, even in the case of the rigid body rotation, the influence of the Coriolis force 
in the reference frame, rotating together with a star, results in the additional frequency splitting which remains in 
the 1-D problem as well. The aim of this work is to evaluate such an effect in a magneto-active plasma beyond the 
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scope of the perturbation theory, apphed in the 3-D model, or JWKB approximation. However we will use the exact 
analytical solutions of the 1-D model. 

This paper is organized as follows. In Sec. [ill we formulate the MHD model for an ideal plasma, discussing all the 
approximations used in Sec. Ill Al In Sec. Ill Bi we evaluate the importance of the differential rotation in RZ. Then the 
approximation of the rigid body rotation is used in RZ. In Sec. lIIII we linearize the complete set of the MHD equations 
and derive the master equation for the z-component of the velocity perturbation Vz{z). With help of this velocity 
component one can receive all the remaining MHD perturbed quantities including the magnetic field perturbations. 
It is the component which results in the Doppler shift of optical frequencies, measured in a helioseimic experiment. 
In Sec, nil Al we consider the case without the magnetic field, accounting for the solar rotation only. 

In Sec, nil B] we derive the simplified master equation for v^^^z) in the limit of relatively weak magnetic field that 
allowed one to avoid the MHD (cusp or Alfven) resonances in the RZ of the Sun. Note that such a resonance was 
examined Refs. [1, Using the magnetic field perturbation theory in Sec. HUB I we obtain the exact analytical 

g- modes spectrum in presence of magnetic fields and solar rotation in RZ. 

In Sec.[lV]we discuss our results. 



II. BASIC IDEAL MHD EQUATIONS 



We start from the following set of MHD equations. The mass conservation law for the total density is of the form: 

(1) 



| + V.(pv) = 0, 



where the total density is the sum of two terms, background density and perturbation, p = pq + p' . 
In the system rotating with the Sun (vq = 0) the momentum conservation takes the form: 



— + (vV)v + 2(0 X v) = 



1. 



P 
-r X 



IP 
8^ 



Sirp 



(HV)H + g - X (rj X r) 



'dt 



(vv)n 



(2) 



where the total magnetic field H = Bq + B' includes the background magnetic field Bq and the magnetic field 
perturbation B', g is the gravity acceleration, Jl is the angular rotation frequency. 
In ideal MHD the Faraday equation takes the form 



'dt 



(H- V)v- (v V)H-Hm, 



(3) 



where m = V • v is the compressibility of gas. 

Finally the full nonlinear ideal MHD equation system is completed by the equation of the entropy conservation: 



d f p 



dt y-p* 



P 



0, 



(4) 



where p — po + p' is the total gas pressure; 7 = Cp/cy is the heat capacity ratio (7 = 5/3 for hydrogen plasma). 
Note, that in the equilibrium condition, 

1 / \ 
V Uo + ^ + g - X (fl X r) = 0, 

the contribution of the centrifugal force (last term) is negligible in the Sun (~ 2 x 10"^). 



(5) 



A. ID rectangular model assumptions 

In order to find spectra of seismic waves accounting for the magnetic field and rotation in RZ a number of as- 
sumptions is required. We shall specify these as they are used, in order to keep clear which results rely on which 
assumptions, but we also list them all here for ease of reference. 
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1. We consider ideal MHD neglecting both the heat conductivity and viscosity contributions to energy losses, as 
well as the ohmic dissipation. 

2. We linearize the MHD equations about an equilibrium background configuration, i.e. a background configuration 
which is time independent and given by the condition ([5]). 

3. We assume the fluctuations to be adiabatic, with the contributions of fluctuations to the heat source vanishing: 

g' = o. 

4. Moreover, we consider the fully ionized ideal gas, that means the thermodynamic quantity, first adiabatic 
exponent 7 = Cp/cy, is time independent and uniform. For numerical estimates we will take 7 — 5/3 for 
hydrogen plasma. 

5. We adopt the Cowling approximation, which amounts to the neglect of perturbations of the gravitational 
potential, (i.e.: (j)' = 0). 

6. We assume the rigid rotation in RZ. 

7. We use a rectangular geometry with the Cartesian coordinates: x, y, and z, where z corresponds to the solar 
radial direction. See Fig. [TJ The background quantities vary along the z direction only (which implies the local 
gravitational acceleration, g, is directed along the z axis, but in opposite direction). We also take a constant 
uniform background magnetic field, Bq, pointing along the x-axis. The rotation (the angular velocity vector $7) 
is assumed to be directed along the y-axis that is perpendicular to eqliptic (xOz-plane) . Far enough from the 
solar center such a slicy ID configuration corresponds to a large scale toroidal magnetic field in RZ in the 3D 
geometry. 



Ay 




z 



FIG. 1: The rectangular (ID) geometry. 

8. The background mass-density profile is assumed to be exponential, po = Pc exp[— for constant Pc and 
H. The conditions of hydrostatic equilibrium for the background then determine the profiles of thermodynamic 
quantities, and in particular imply 7 is a constant. 

9. We assume that Brunt- Vaisala frequency is zero in CZ and non-zero, but constant, in RZ. 

This last assumption applies to very good approximation for real mass-density profiles obtained by Standard Solar 
Models, provided we identify the z direction with the radial direction, and focus our attention to deep within the 
radiative zone. The constancy of 7 in this region is also expected since the highly-ionized plasma satisfies an ideal gas 
equation of state to good approximation. The rectangular geometry provides a reasonable description so long as we 
do not examine too close to the solar centre. What is important about our choice for Bq is that it is slowly varying 
in the region of interest, and it is perpendicular to both g and all background gradients, Vpo, Vpo, etc. 

As suggested in Q such ID picture can be fully described in analytical terms in contrast to the 3D case. 

B. Differential rotation 

In 3D-geometry (spherical coordinates), for which = cos 6*6,. — sin 9eg, in the absence of angular acceleration, 
dVt/dt = 0, the last term in Eq. ^ , — rsin0(vV)rie$, is much less than the Coriolis term in the l.h.s. of Eq. ^ 
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under the condition 



or 



r sin I Vr + — 1 *C 2f2 [sin 0Vr + cos 6'u6il . 
^ dr r 80 ' 



dfl dfl 
r— <2r2, tan6i— <2r2. (6) 
ar at' 

There are no doubts that such conditions can be fulfiUed in RZ (for almost rigid rotation) while in CZ the last 
inequality is valid near equator 9 = n/2, (inclined on 7 degrees with respect to the eqliptic latitude d — 0°), only in 
the limit 9f2/96'|e^^/2 0- Looking at the solar rotation profile in Fig. 1 (right panel) in Ref. fu\ demonstrating 
the SoHO results we see that the first condition in ([6]) is valid at least for low and middle latitudes with an accuracy 
- An/n - 50/400 = 0.125 < 1. 

In our rectangular geometry the corresponding conditions: 

dfi „^ dfi ^^Vx 
z— <217, x—^2n—, 7 

az az Vz 

are similar to ^ where the velocity ratio is larger than the unity, v^/v^ > 1, especially at the MHD resonances where 

In what follows we put = const, or neglect angular acceleration and the differential rotation omitting the last 
line in Eq. 

III. LINEAR IDEAL MHD IN ROTATING SUN 

Assuming the plane wave propagation along the background magnetic field, ^ ^-ti^t-k^x) ^ linear ideal MHD we 
obtain from the full MHD system (jlj- ([4]) the following set of eight (=8) equations: 

J _ Vz dpo ^ po ( dvy ^ dvz\ , poVx 



= + + ^ (8) 

ujj (iz luj \ ay az J V 



dP . , ikxBo 



ipoLoVy + ^^S;, (10) 



— = ipoWVz H Bz + 2pQilVx - pg, (11) 

oz 47r 



J/ _ Bo f dvy ^ dvz 



B'x^-{Y + ^ I (12) 
luj \ oy az J 



K = (13) 



B'z = -^vz, (14) 



V 



, 1 

ci 



p ^ BqBx 

47r 



Vz (7 - 1) dpo 
iuj 7 dz 



(15) 
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where V = uj/kx is the phase velocity of MHD oscillations we are looking for; P = p' + BqB'^/Att is the total pressure 
perturbation in linear approximation; Cg = \/'jgH is the sound velocity; H — O.IRq is the density scale height in the 



background density po ~ PcG 



-z/H 



Pc 



150 g/cm is the central mass density in the Sun. 



Using Eq. 
components: 



Q, from which one finds P — BqB'^/Att, then equating ([5]) and ([TS]) we derive in terms of 



1 / dVy 

dy 



dvz 
dz 



dpo 2Q.Vz 



iuj^Po dz ikxC^g 



(16) 



We see from Eq. (jlOl) that in our geometry we can put Vy — By — d /dy — 0, or RZ medium is uniform along rotation 
axis (y-axis). This allows us with the use of Eq. ^ for the pressure P, Eq. (|12|) for B'^ and from Eq. (fT6|) to 
obtain from Eq. pT|) the master differential equation of the second order: 



dC 



i + r — 



1,2 „2 



2Hn 



2 

— 1 

,7 



(17) 



Here we input the Brunt- Vaisala frequency N , N"^ = {cl/ H'^){j — l)/7^ = {9/H){-y — l)/7, on which we can normalize 
the compressibility parameter ai = uj'^/kxcl = .2A{uj / N)"^ / where the coefRcent 0.24 = (7 — l)/7^ is given by the 
heat capacity ratio (7 = 5/8) for hydrogen plasma. The argument C = (1— £, = Coe^^'^, is given by the background 
magnetic field Bq, = k'^v'j^Q/uj^, where vao is the Alfven velocity at the center of the Sun, vao — (-So/43.4 G) cm/sec, 
given by the central density pc- 

The solution of this equation can be easily obtained from the Gauss equation for Vz — C'^^(C) i^^ analogy with [sj: 



Vz{(:) = DirFia,b;c;C)+D2C 



where F{a, b; c; C) are Gauss (hypergeometric) functions, 



F{1 



l + 6-c;2-c;C), 



(18) 



(^1,2 




- 1 



4f}2 



217/2 ^ 

utk \7 



(19) 



are the upperscripts in the solution (upper sign for ai), a — ai + K , b — ai ~ K , c — a + b, K — k^H. 

Using the resonable boundary conditions (see e.g. in [8|) one can derive the MHD dispersion equation for the 
magneto-gravity modes u!{K, Bo,n,) expressed through hypergeometric functions in at the fixed point (at the 
center of the Sun, z — 0, ( ~ = {1 — ai)^Q). However, the analysis of such dispersion equation is still very 
complicated. Below we consider some simplifications of the master equation p?)) which allow us to get analytic MHD 
spectra in the RZ. 



A. Zero magnetic field limit 

In the limit C = (1 — ai)^ 0, Bq ^ 0, with the use of the transformation Vz = e^^^^'i'^z) and the change 
Cd/dC = Hd/dz, Eq. ^ converts to: 



d^ 

dz^ 



Am 



2 

1 

7 



2n 

ujk^H 



(20) 



In the absence of rotation, = 0, this equation is similar to 3D osillation Eq. (7.90) in (T3 |: 



d^ 

dr^ 



1 



4i72 



- 1 



^ = 0. 



(21) 



Rewriting Eq. ^ for f7 = as 



^ + ^*(z) = 



dz2 
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where /?o = ^/IlO{NyII^^^^Ty~n'TAlO^ in RZ, and as 



in CZ, where Ncz — 0, or Pq ^ iTo = i^J AK'^{1 — oi) + 1, then matching the corresponding solutions at the top of 
RZ, z = zjiz, one can easily obtain the g-mode spectrum in our ID model 



2H 



Tin = arccot 



(22) 



Here n = —1, —2, —3, .... The analogous solution of Eq. (PO)) for the rotating Sun, il ^ Q, leads to the spectrum: 

Pnznz 



2H 



Tin — arccot 



'f3n 



where 



\ 



0.96 



2n 



1 + 4fc2 



iV2 



1 /2r2 



5k \ N 



(23) 



(24) 



To 



/4fc2 



;2 - 4f^2 



1 
5^ 



20 



(25) 



Here n = —1,-2,,... (see Fig. [5]). We can neglect a term in ^2 Eqs. ([M]) . (|25p . since the rotation parameter 
2n/N = 2 X 10^'^ is small comparing with the frequencies in the GOLF experiment 1 > uj/N > 0.2, we are 
interesting in. 

The rough estimate of the rotation frequency splitting w — s- w — Auj due to the Coriolis force near lu < N, can be 
easily found from Eq. accounting for the shift of the parameter N'^ /oj'^ [N"^ /uo'^){\ + 2/S.u)/lo): 



an 



UJ 



1 



IQK 



20 \ 



)(^) 



1 

5A' 



X 10" 



(26) 



Thus, for isT > 1 and waves Vz propagating in eqhptic (where g _L fJ) this frequency shift is negligible, an ~ 10"''. 

In Fig. [2] we show the g-mode spectrum Lj{n^K) accounting for rotation O in the Sun and given by the dispersion 
equation (j23p for a few modes, n — —1, —5, —10. In Fig[3]we show g-vaoAe frequency shifts ai^{K) — (lj — a;o)/^o in 
dependence of the transversal wave number K neglecting RZ magnetic field, Bq = 0. Here the frequency responds 
to the absence of rotation, O = = 0. One can see that the rotation correction originated by the Coriolis force 
increases if K decreases, especially for the highest frequency n = —1. This behaviour is in an agreement with what 
one expects in the 3D case (see discussion below). 



B. Magnetic and rotation corrections to the g-mode spectrum 



In order to obtain spectra of g-modes in the presence of the magnetic field let us define the coefficent in front of 
the second derivative in Eq. ([T7]) as 1 — C, where C, = v\{l — ai)/V^. 

We consider the perturbative regime for magnetic fields, where Alfven velocity is much less than phase velocitiy, 
v\ = v^^e""/" ^.V"^ = Lo'^/kl = [N"^ H"^ I K'^]{u} / , so that C < 1 and MHD resonances (( = 1) do not appear 
within the RZ (see Fig. 2 in 5]). The perturbative magnetic field for which the maximum Alfven velocity is small, 
v\{Q.7Rq) < y2^ obeys 

«„«|(^)Ma 



This condition constrains the range of magnetic field strength when the perturbative approach is valid. One can check 
that this condition is well satisfied for the magnetic fields and for the modes down to n = —10 considered here. 
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FIG. 2: The g-mode spectra accounting for the solar rotation versus wave number, K, shown for different mode numbers, 



,x 10 




FIG. 3: The relative shift of g-mode frequency due to rotation, {u — uoq)/uo, versus wave number, K, shown for different mode 
numbers, n — —1, —5, —10. 



The absence of MHD resonances allows us to set 



-1/2 



and thus derive from the master Eq. (fT7|) the 0(?;^/F^)-correction to Eq. ([201 



d^ 



PI 



* = 0, 



where we defined f3n in Eq. (|24|) and the wave number correction kj^{fl) is given by 



(m 



0.24 



Introducing the coefficent 



Ao - (1-ai) 



T/2 



Ar2 



mV- /0.24 



1 /2r2 
5^ 



1 /2ri\ /TV 



(27) 



(28) 



(29) 
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and using the change of the variable 2s 



z/H + \n{4:K^Ao 



[A 



), one can rewrite Eq. ([27)1 as 
2^1 *(s) = 0. 



(30) 



This equation has general solution expressed in terms of Bessel functions of the first kind, ^'(s) = CiJiii^{e^) 



C2J-i[3a{^^) where Ci,2 are constants, 
boundary condition at the solar center, 
equation ([3D]) in RZ, < z < zrz'- 



Then following Q we can find the corresponding solution in RZ using the 
^'(z = 0) = (as Vz{0) = 0), excluding one constant to obtain the solution of 



J.,p^{2KAy^) 



J. 



i^KA^o 



l/2z/2H 



(31) 



In CZ the Brunt- Vaisala frequency N vanishes, iV = 0, so that we are led to the same solution as in isotropic case, 
^cz{z) — Ccz exp(— zF^/i/), where is given by Eq. (^5]) . and we safely neglected the possible influence of a CZ 
magnetic field (see comment on that after Eq. (25) in 5]). Since the arguments of Bessel functions are small, we can 
use the first terms in the Bessel function series only. 

Then by matching the logarithmic derivatives of the solutions ^ rz and ^cz at the top of RZ, z — zrz, one obtains 
the dispersion equation for the case _Bo 7^ which generalizes our result in Eq. (27) f5l| accounting for rigid rotation 
of RZ in the Sun: 



1 



X cot 



1 + A 

zrz 
H 



I 



1 + 0, 



- Pf^ sm{l3nZRz I H) 

-Ff2, 



(32) 



where krz — KA^'^e^'^^^^^, cocfficents Aq, Pn and Fq are given by Eqs. ((29)) . pi)) and ([25)) respectively. 

This is our main equation giving the g-mode spectrum wb{K) in ID MHD model accounting for the rotation of 
the Sun. In the Fig. [4] we illustrate the rotation splitting originated by the Coriolis force for the magnetic field 
Bo = 700 kG in dependence on the transversal wave number K for a few modes n = —1, —5, —10 (dashed lines there). 



0.012 




1.5 1 1.5 2 2.5 3 3.5 4 
K 



FIG. 4: The relative shift of g-mode frequency due to rotation and magnetic field, versus wave number, k, shown for magnetic 
field Bo = 700 kG and diflPerent mode numbers, n = —1, —5, —10. Solid lines - neglecting rotation, dashed lines - accounting 
for the rotation. 



IV. DISCUSSION 



As we expected the contribution of the Coriolis force to the g-mode frequency shift is small: even for the radial 
order n = — 1 one obtains ao < 0.003. In Fig. [5] we show how the magnetic field shift as [l] is corrected due to the 
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presence of the solar rotation, 7^ 0. We compare there two curves calculated for the magnetic field Bq = 700 kG: (i) 
as (if) = {cobo — '^o)/'^o given by the solid line and corresponding to the previous result obtained in in the absence 
of rotation, = 0, and (ii) aQ{K) — {lub ^ ^o)/^o given by the dashed line where ujb accounts for the magnetic field 
and the Coriolis force contributions. 

Note that in our ID approach, in which azimuthal number m is absent and the transversal wave number K is 
analogous to the angular degree I (the Legendre number), the decrease of this wave number leads to an enhancement 
of the Coriolis force contribution, as it can be seen on Fig. [H in the rotational splitting as it is predicted in the 3D 
case. 

Let us note also that a crude approximation for the Brunt- Vaisala frequency profile N{z) in RZ relies here on 
a fixed value A'' = constant for which we have chosen the value = 10"'^ rad/sec as in Ref. 0. This value was 
used in all calculations instead of the maximum value of the Brunt- Vaisala frequency, N = 2.8 x 10^'^ rad/sec 
{vn ~ N/2tt ~ 440 /iHz in cycling frequencies, see, e.g., in Ref. [i2|)- Relying on such maximum N one gets 
2.8^ — 7.8 times smaller effect for the same value of magnetic field. And, vice versa, the same magnetic shift of 
the order 1 % (together with the Coriolis force correction) can be obtained relying on the Brunt- Vaisala frequency, 
= 2.8 X 10~^ rad/sec, if a stronger (2.8 times greater) magnetic field exists in RZ. For example, to receive the 
corresponding frequency shift for the g-mode with n = — 10 the magnetic field should be Bq — 700 kC 

Therefore magnetic fields resulting in the shift Su/u! ~ 1% appear to be quite strong, in particular for high frequency 
g-modes, close to the Brunt- Vaisala frequency, observed in the experiment. Thus for the corresponding shift of the 
mode with 1/ — 222.46 /iHz, one predicts the strong RZ magnetic field B > 8MG (see Ref. (gJ) that already exceeds 
the existing central magnetic field limits. As it follows from the solar poles oblateness ~ 10^^, one gets the limit 
B < 8MG (see Ref. [l^). This discrepancy can be explained by the fact that the observed mode shift results from 
its interaction with other g-modes, close to the Brunt- Vaisala frequency, rather than a magnetic field or a rotation. 

The configurations with 3- vectors g, $7 and Bq, different from that shown on Fig. [1] in the same 1-D rectangular 
geometry could lead to the different spectra we do not discuss here. 

In conclusion we mention that the ultimate goal of all solar studies would be the complete 3D MHD description of the 
solar evolution. However, it is rather far future from now. On the way to this goal in this letter we present qualitatve 
analytical ID model which describes the influence of RZ magnetic field on helioseismic wave spectra accounting for 
the solar rotation. 
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